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Abstract 

We investigate the use of partially twisted boundary conditions in a lattice 
simulation with two degenerate flavours of improved Wilson sea quarks. The 
use of twisted boundary conditions on a cubic volume (L 3 ) gives access to 
components of hadronic momenta other than integer multiples of 2n/L. Par- 
tial twisting avoids the need for new gluon configurations for every choice of 
momentum, while, as recently demonstrated, keeping the finite-volume errors 
exponentially small for the physical quantities investigated in this letter. In 
this study we focus on the spectrum of pseudo scalar and vector mesons, on 
their leptonic decay constants and on Zp, the matrix element of the pseudo 
scalar density between the pseudo scalar meson and the vacuum. The results 
confirm the momentum shift imposed by these boundary conditions and in 
addition demonstrate that they do not introduce any appreciable noise. We 
therefore advocate the use of partially twisted boundary conditions in appli- 
cations where good momentum resolution is necessary. 

PACS: 11.15 Ha, 12.38 Gc 



1 Introduction 

In lattice simulations of QCD on a cubic volume (V = L 3 ) with periodic boundary con- 
ditions on the fields, the components of hadronic momenta pi are quantized in integer 
multiples of 2n/L. For currently available lattices this implies that the lowest non-zero 
momentum is large, typically 500 MeV or so, and there are large gaps between neighbouring 
momenta. This limits the phenomenological reach of simulations, particularly for momen- 
tum dependent quantities such as the form-factors of weak semileptonic decays of hadrons. 
In ref. [1] Bedaque proposed the use of twisted boundary conditions 1 for the quark fields ip 

ip(xi + L) — e i6i i/j(xi) . (1) 

Twisted boundary conditions allow for simulations with arbitrary components of hadronic 
momenta. For example, the momentum of a meson composed of a quark with flavour 1 

— # 

satisfying boundary conditions with a twisting angle 9i = (9 n , 6* 12 , #13) and an antiquark 

— # 

of flavour 2 with angle 9 2 is 




where n is a vector of integers. 

The practical difficulty in using twisted boundary conditions in lattice simulations with 
dynamical quarks is that it requires the generation of a new set of gauge field configurations 
for every choice of twisting angle(s). In refs. [3,4] it was shown that for many physical 
quantities one can use partially twisted boundary conditions, i.e. impose twisted boundary 
conditions for the valence quarks but periodic boundary conditions for the sea quarks, 
thus eliminating the need for new simulations for every choice of momentum and making 
the technique practicable. The physical quantities for which partially twisted boundary 
conditions can be applied include those with at most a single hadron in the initial and final 
states (and possibly even in intermediate states), for which the finite- volume effects de- 
crease exponentially with the volume. For these processes the finite-volume effects depend 
on the twisting angle(s) but remain exponentially small. 

For some processes with energies above a two-body threshold, such as K — > nn decays with 
the two-pions in an isospin zero state, the finite-volume effects decrease only as powers of 
the volume and must be subtracted for acceptable precision to be reached. We are not 
able to perform these subtractions if partially twisted boundary conditions are used. Here 
we will only consider processes for which such a problem does not arise. 

In this letter we confirm the theoretical results of ref. [3] in a numerical study of partially 
twisted boundary conditions for dynamical, non-perturbatively improved Wilson fermions. 
In particular we find that: 

• The energies of n and p-mesons (with masses below the two-pion threshold) satisfy 
1 See the references cited in [1] and [2] for earlier related ideas. 
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the expected dispersion relation 
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where Q\ and 62 are the twisting angles of the two valence quarks and p\ at = (2ir/L)n 
is the contribution to the meson's momentum introduced by the Fourier transform 
of the correlation function. 

This study extends the one in ref. [2] where the dispersion relation for pseudo scalar 
mesons with twisted boundary conditions in the quenched approximation was found 
to be consistent with expectations. 

• The values of the leptonic decay constants of 7r and p mesons and of the matrix 
element (0|P|7r) of the pseudo scalar density P are independent of the twisting angles 
as expected. 

A further reassuring result of our study is that twisted boundary conditions do not in- 
troduce additional noise in the data. As we increase the meson's momentum by suitably 
varying the angles 6\^, the statistical errors on meson masses and matrix elements increase 
smoothly. However, when comparing results obtained with twisted and periodic boundary 
conditions with similar momenta (i.e. momenta close to 2n/L or y / 2(2vr/L)) the errors are 
found to be comparable. 

The plan of the remainder of this letter is as follows. In the next section we present 
the details of our computation, the parameters of the simulation (including the choice of 
twisting angles) and a description of the analysis. We present our results in sec. El and 
conclusions in sec. EJ 



We study meson observables on sets of gauge configurations which were generated with 
two degenerate flavours of sea quarks using non-perturbatively improved Wilson fermions 
and the plaquette gauge action on the torus with periodic boundary conditions (/3 = 5.2, 
a « 0.1 fm, c sw = 2.0171, {L/af x T/a = 16 3 x 32). We used the ensembles of field 
configurations which were studied in detail in [5, 6] and took over the suggested separa- 
tion of measurements by 40 trajectories in our analysis. The simulated quark masses are 
summarized in tabled Propagators and correlators were calculated using the FermiQCD 
libraries [7-9]. We stress that the aim of the present study is to investigate the consistency 
and effectiveness of using partially twisted boundary conditions at fixed values of the quark 
mass. We do not attempt to perform a chiral extrapolation. 

For each flavour of valence quark we impose the boundary conditions in eq. for a variety 
of twisting angles 6 = {61,62, 6 3 ). When evaluating the corresponding propagators we make 
use of the change of quark field variables 



2 Details of the Simulation and Analysis 
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N 

1 ' mcas 


0.13500 


0.697(11) 


200 


0.13550 


0.566(16) 


200 



Table 1: Simulation parameters. 



where ip(x) satisfies periodic boundary conditions. The phase factor cancels in all terms 
of the lattice fermion action except for the spatial hopping terms which now become (for 
2 = 1,2,3) 

4>(x) [e^K(x)(l - 7i )^(x + + e'^Ujix - i)(l + n)$(x - i)] . (5) 

In practice therefore, the partially twisted quark propagator can be computed by inverting 
the standard improved Wilson-Dirac operator in a gauge field background where the link 
variables {Ui(x)} have been replaced by {e l ~^Ui(x)}. 

The physical observables which we study in this letter are the energies and leptonic decay 
constants of the pseudo scalar and vector mesons and the matrix element of the pseudo 
scalar density. In order to determine these, we compute the following correlation functions: 



C Ao p(t,p) 


= ^e^-*<0|4(£, 0^(0)|0>, 

X 

= ^e^^OlPOr^P^lO), 

X 


(6) 


C PP (t,p) 


(7) 


C Ao A (t,p) 


= ^e^->|4(f,t)(4(0))t|0), 

X 


(8) 


C ViVi (t,p) 


= ^e^ lat ^(0|y/(f,t)(y/(0)) t |0) (no sum on i) , 


(9) 



where P(x) is the pseudo scalar density 

P(x)=^ 2 (x) l5 Mx) (10) 

for quarks of flavour 1 and 2 (with twisting angles 6\ and 62), and Vj(x) and A T (x) are 
the improved vector and axial- vector currents 

v i( x ) = ^i{ x )l^i{ x ) + ac v (g ) hydl + d v )ip 2 (x)a IJU ,ip 1 (x) 
A K X ) = ^2( x h^5M x ) + ac 49o)^(d* + d^)P(x). 

Here, and d* are the forward and backward derivatives and cv(go) and ca(<?o) are im- 
provement coefficients which we take from [10] and [11] respectively. Since we are primarily 
interested in the effects of twisted boundary conditions we do not attempt to compute the 
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renormalization constants of P, and A^, nor do we implement improvement factors of 
the form 1 + b(go)m q a, where m q is the mass of the quark. The inclusion of these factors 
would of course be necessary if we were attempting to determine the physical leptonic 
decay constants. However, they are overall factors for each choice of quark mass and are 
independent of the twisting angles, while it is precisely the dependence on these angles 
which is the object of our study. 

The momentum, p, of the meson is given by 

P = Plat T , (11) 

where p\ at = (27r/L)n and n is a vector of integers. 

At large values of t the time dependences of (jOJ-Q approach: 

C Ao p(t,p) - -^rZ P M (p)e-^ T / 2 8mh((t-T/2)E n ), (12) 

C PP (t,p) -> ^Z|e-^ r / 2 cosh((t-T/2)K), (13) 

C AoAo (t,p) -> yj- M 2 (p) e~ E * T l 2 cosh((t - T/2)E n ), (14) 

C ViVt (t,p) - -^N?(p)e- E ^ 2 cosh((t-T/2)E p ) (< = 1,2,3), (15) 

where, for each choice of quark masses, we have denoted the lightest pseudo scalar and 
vector mesons by 7r and p respectively and and E p are the corresponding energies which 
we expect to satisfy the dispersion relations in eq. 0- The notation for the matrix elements 
is as follows: 

Z P = (0|P(0)|tt(p)>, (16) 
M (p) = (0\A (0)\7r(p)} = f 7T E 7T} (17) 

nKp) = Ei< i^( )i^ A ))i 2 = /x( 1 + ^) ( 18 ) 

where the index A labels the p-meson's polarization state. 

In this letter we study the validity of the dispersion relation in eq. and the independence 
of f n , f p and Z P of the momentum. We evaluate the quark propagators for four values of 
the twisting angle 9: 

9 = 6, (2, 0, 0), (0, 7T, 0) and (3, 3, 3) . (19) 

For each value of /t va i, quark and antiquark propagators with all possible pairs 9\ and 9i 
were combined to construct correlation functions for mesons with a variety of momenta. 
Moreover we also combined them with Fourier momenta pi at = (0, ±27r/L, 0) to increase 
the range of momenta which can be reached. When presenting our results in the following 
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section, we include for comparison results without twisting (0 X = 2 = 0), obtained by 
averaging over the 12 equivalent momenta with |pi at | = a/2 x 2tt/L and those obtained 
by averaging over the eight equivalent momenta with |pi at | = a/3 x 2tt/L. Of course this 
averaging reduces the statistical errors and this should be borne in mind when comparing 
the errors at these untwisted momenta with those at momenta with 9\ — 62 / for which 
such averaging is not possible. 

Applying the jackknife procedure to the data for the correlation functions in eqs. 
we have extracted all observables in the pseudo scalar channel from a combined non-linear 
X 2 fit to the functional form suggested by (JUJ-fHU), (fTtjj) and (fTTj) . The fit-ranges were 
chosen to yield compatible results under variation of the range by at least one unit in 
t/a. We applied the same procedure in the vector channel, combining the data for the 
correlation function for i = 1,2 and 3 in one fit using the expressions in eqs. (fTH|) and 

3 Results 

The series of plots in figures and El show our data as a function of (pL) 2 in the range 
\pL\ G [0, a/3 x 27r]. Fig. ^contains the results for the energies as a function of momentum 
and fig. El those for the decay constants and Z P . To ease orientation, the positions of the 
discrete Fourier momenta |pi a t-k| = 0, 27r, a/2 x 27r and a/3 x 27r are indicated by dashed 
vertical lines. We emphasize that it is only at these values of momenta that one can obtain 
results using periodic boundary conditions. In fig. [2] we zoom into the region \pL\ < 2ir for 
the dispersion relations. In this region we would expect lattice artefacts to be small and 
the use of twisted boundary conditions to be particularly useful. 

In each plot, the (blue) triangles correspond to points in which the correlation function was 
evaluated with p lat = 0, but with all possible pairs of 9i and 9 2 from the set in (JTHJ). The 
(red) diamonds and (green) squares represent the results obtained with pi at = (0,27r/L,0) 
and pi a t = — (0, 2tt/L, 0) respectively, combined with all possible pairs of d\ and 62- The 
four points with 9\ = 62 = with \p\ a t\ = 0, 2tt/L, a/2 x 27r/L and a/3 x 27r/L are denoted 
by (black) circles. 

For the discussion of our results it is convenient to rewrite the dispersion relation in eq. (JHJ) 
in the form 

(aE n/p ) 2 = (am Vp ) 2 + A 2 (pL) 2 (20) 

where A 2 = (a/L) 2 = 0.0039. The dispersion relation (J2"U|) is displayed as the dashed line 
in the plots of fig. In the first row of table |2] we present the xVd.o.f of the comparison 
of our data to eq. ()20)1 over the range < \p \ 2 L 2 < (27r) 2 using the values of the meson 
masses obtained from fits at zero momenta. In the third row of the table we present the 
values of A 2 obtained by fitting the lattice data to the functional form in eq. (J2*U|) over 
the same range in momentum, but allowing A 2 to be a parameter of the fit. We note that 
our values for the ratios m n /m p agree with those found earlier on the same configurations 
in [5]. 
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Figure 1: The plots in the first line illustrate the results for the dispersion relation for the 
7r and the p (empty and full symbols respectively) for the two choices of the quark mass. In 
the second line we show the corresponding relative error as a function of the momentum. 



0.13500 k = 0.13550 



7T 


P 


7T 


P 


0.3 


1.0 


0.6 


1.7 


1.8 


2.5 


0.9 


2.1 



X 2 /d-o.f|j23 

rVd.o.figxj 

A 2 from H2) 0.0040(1) 0.0042(2) 0.0040(1) 0.0048(4) 

Table 2: x 2 /d.o.f. for the lattice data with respect the expectations eqs. (|2T)|) and ()21)1 with 
A 2 = a 2 /L 2 = 0.0039 (first two rows) and the results obtained from a fit to (J20J) with A 2 
left parameter of the fit (third row). 
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Figure 3: The first line shows the results for the ir and p decay constant (empty and 
full symbols respectively) and the second line shows the matrix element (fTK|) for the two 
choices of the quark mass. In each plot the horizontal lines represent the central value at 

Plat = £ =^2 = 0. 
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As the momentum of the meson grows so do the expected discretization effects in the 
dispersion relation. For example, a free scalar particle with a wave-function 4>(x) satisfying 
the (Minkowski-Space) Klein-Gordon equation (□+m 2 )0(x) = 0, with a generic discretized 
second derivative defined by d 2 f (x) / dx 2 = (f(x + a) + f(x — a) — 2f(x))/a 2 , satisfies the 
following lattice dispersion relation: 

srf (^) =sinh2 (^) +Esin2 (A^). ( 21) 

i 

To indicate the possible size of lattice artefacts we plot the dispersion relation of eq. (|2*T|) as 
the solid curve in fig. ^ We stress, however, that in an interacting theory the discretization 
errors will in general be different from those in eq. (|21|). Indeed there seems to be no 
evidence from our data that eq. (|2*T|) is a particularly good representation of the lattice 
artefacts (see the second row of table El where we show the x 2 /d.o.f. from a comparison of 
our data with (|21jl ). At small momenta the solid curve merges of course with the dashed 
line representing the continuum dispersion relation. 

We conclude from the results for the dispersion relation plotted in figs. [U and El that the 
use of partially twisted boundary conditions is beautifully consistent with expectations, 
particularly at low momenta where the lattice artefacts are small 2 . An important further 
observation is that there is no evidence in our data that the introduction of twisted bound- 
ary conditions increases significantly the statistical or systematic uncertainties. This is 
illustrated in the second line of figure which shows the relative error in the pion energy 

Se. = ^ (22) 

as \pL\ is varied. The error SE n is the jackknife error, including the statistical error and the 
systematic uncertainty stemming from the improvement constants in the improved quark 
currents. The plot shows that the errors increase smoothly as the momentum increases, 
and that no appreciable additional noise is introduced by partial twisting 3 . We observe 
the same behaviour for all analyzed quantities. 

In fig. El we plot our results for the decay constants f n and f p and for Z P . The values 
for af n agree with the ones obtained in [6] at \pL\ = 0. Again we see that the results 
are completely consistent with theoretical expectations, being independent of the twisting 
angles and Fourier momenta. 



2 The different data points at \pL\ = and 2tt correspond to 0\ — 62 but with different choices of d\ 
and 82 ■ 

3 For k = 0.13550 at \pL\ = 2ir/L we observe a fluctuation in the effective mass from one of the gauge 
configurations and with our choice of the position of the source (this fluctuation has also been observed by 
our colleagues in the UKQCD collaboration [12]). The fluctuation is particularly noticeable when twisting 
both quarks by 9 — (3, 3, 3) and this is the reason for the larger jackknife error at this particular momentum 
and twist combination (see figure 
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4 Conclusions 



We have investigated the use of partially twisted boundary conditions in evaluating the 
energies of pseudo scalar and vector mesons and their leptonic decay constants. The 
results are very encouraging; it does appear that the method allows the evaluation of 
physical quantities with any momentum. Moreover the use of these boundary conditions 
does not appear to increase the errors in any appreciable way. It will be important to 
monitor whether this continues to be true as the quark masses are decreased. Once the 
quark masses are such that two-pion intermediate states contribute significantly to the p- 
meson's correlation function the finite-volume effects will no longer fall exponentially with 
the volume, but only as powers. For the pion observables studied in this paper this is not 
the case. 

Partially twisted boundary conditions will be particularly useful for evaluating momentum- 
dependent physical quantities. One important application is to the determination of the 
form-factors of semileptonic weak decays of heavy (D and B) mesons to light mesons. 
For these processes, with conventional periodic boundary conditions, the initial and final 
state hadrons are restricted to have momenta (2n /L)n where n is a vector of integers. In 
order to avoid lattice artefacts the possible values of \n\ are frequently limited to 0, 1, 
and perhaps Thus, for any particular choice of quark masses, the number of values of 
the momentum transfer, q 2 , or the light meson energy, E, is also very limited. Moreover, 
chiral extrapolations are conveniently performed at fixed q 2 [13-15] or fixed E [16-19] (and 
heavy quark extrapolations at fixed E), but q 2 and E vary with both the momentum and 
quark masses. Ansatze for the form factors, such as the Becirevic-Kaidalov [20] model, 
are used to interpolate and extrapolate simulation data to sets of common q 2 or E values 
before the extrapolations are performed. Using twisted boundary conditions would enable 
the form factors to be evaluated directly at these common values, removing the need for 
the intermediate form-factor fit. 

For some other physical quantities, such as the moments of hadronic deep inelastic structure 
functions or light-cone distribution amplitudes, it may be helpful to use twisted boundary 
conditions even though it is not strictly necessary. The corresponding matrix elements are 
proportional to factors of pi, where p is the momentum of the hadron, so that the correlation 
functions must be computed with p ^ 0. The use of twisted boundary conditions allows 
\p | to be decreased and hence the lattice artefacts to be reduced. Moreover by varying p 
one can verify that the leading twist component has been extracted correctly. 

Following the successful conclusion of this exploratory numerical study of the implemen- 
tation of partially twisted boundary conditions we now look forward to applying them in 
lattice computations of a wide variety of phenomenologically important quantities. 
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